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ABSTRACT 

We have obtained a moderately long (10 5 second) ASCA observation of the Seyfert 
1 galaxy Markarian 231, the most luminous of the local ultraluminous infrared galaxy 
(ULIRG) population. In the best-fitting model we do not see the X-ray source directly; 
the spectrum consists of a scattered power-law component and a reflection component, 
both of which have been absorbed by a column Njj ~ 3 x 10 22 cm -2 . About 3/4 of the 
observed hard X-rays arise from the scattered component, reducing the equivalent width 
of the iron Ka line. The implied ratio of 1 — 10 keV X-ray luminosity to bolometric 
luminosity, L x /Lj, i ~ 2%, is typical of Sy 1 galaxies and radio-quiet QSOs of comparable 
bolometric luminosities, and indicates that the bolometric luminosity is dominated by 
the AGN. Our estimate of the X-ray luminosity also moves Mrk 231 in line with the 
correlations found for AGN with extremely strong Fe II emission. A second source 
separated by about 2 arcminutes is also clearly detected, and contributes about 25% of 
the total flux. 

Subject headings: galaxies: active - galaxies: individual (Markarian 231) - galaxies: 
nuclei - galaxies: Seyfert - X-rays: galaxies 
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1. Introduction 

Sixteen years ago, the IRAS all-sky survey revealed the existence of a class of objects which 
emit nearly all of their energy at far-infrared wavelengths (Soifer et al. 1984). At the high- luminosity 
end of the distribution, Ljr > 10 12 L , the so-called "ultraluminous infrared galaxies" (ULIRGs), 
the space density of such objects is approximately twice that of optically-selected QSOs, the only 
objects with comparable bolometric luminosities (Sanders & Mirabel 1996). The far-infrared energy 
distributions of the prototype objects, Arp 220 and NGC 6240, are very similar to that of the classic 
starburst galaxy M82, peaking at A ~ 60/xm, and do not exhibit the "warm" IRAS colors typical of 
Seyfert galaxies (DeGrijp et al. 1992). A large fraction of the ultraluminous galaxies show evidence 
for interaction or merging with a companion galaxy (Sanders et al. 1988; Melnick & Mirabel 1990), 
and most show luminous molecular line emission, indicating substantial quantities of molecular gas, 
which could potentially fuel either a burst of star formation or an AGN (e.g., Sanders & Mirabel 
1996, and references therein). 

Ever since the discovery of this class of galaxies debate has raged over the nature of the 
dominant source of energy. The failure to detect any of the ULIRGs at X-ray wavelengths with 
HEAO-1 (Rieke 1988) strengthened a starburst interpretation, as did the recent failure of the OSSE 
instrument on the GRO to detect Arp 220, Mrk 231, or Mrk 273 in the 50 - 200 keV energy range 
(Dermer et al. 1997); this rules out very luminous Compton-thin sources in these galaxies. Most 
recently, observations of a number of ULIRGs with the short-wavelength (SWS) spectrometer on 
the Infrared Space Observatory failed to detect any high-excitation lines in the mid-infrared, as 
might be expected from even a dust-shrouded AGN (Sturm et al. 1996; Genzel et al. 1998), and 
this has been taken to be evidence for starbursts rather than AGN. However, it is now clear that in 
many ULIRGs the column densities toward the nucleus are so high that the dust optical depth is 
greater than unity even in the mid-IR. Furthermore, the discovery of Compton-thick X-ray sources 
in a number of ULIRGs, notably including NGC 6240 (Iwasawa & Comastri 1998; Vignati et al. 
1999), substantially weakens the argument against AGN-dominated sources based on hard X-ray 
non-detections. Such sources, in which even hard X-rays are seen only through reflection, are 
characterized by a relatively flat continuum above E ~ 3 keV and very large equivalent widths of 
the iron Ka line, exceeding 1 keV (e.g., Krolik & Kallman 1987; Krolik, Madau k, Zycki 1994; 
Matt, Brandt, & Fabian 1996). Such large equivalent widths can be produced only if the X-ray 
continuum incident on the line-emitting region is much larger than that which we observe. If the 
reflected fraction in NGC 6240 is similar to that in NGC 1068, the bolometric luminosity of the 
AGN is of order 10 45 erg s _1 , and dominates the far-infrared output of the galaxy. 

Foremost among the local ULIRGs is Markarian 231 (z = 0.042) , which is the most luminous 
object in the local universe (z < 0.1) (Soifer et al. 1984). Its bolometric luminosity is L^ i ~ 
3 x 10 12 L (H q = 75 km s" 1 Mpc -1 ), nearly all of which emerges in the far-infrared. Although 
classified as a Seyfert 1, Mrk 231 shows evidence for substantial amounts of reddening (Ay ~ 
2: Boksenberg et al. 1977) and when corrected for extinction its optical luminosity approaches 
quasar levels. Furthermore, Mrk 231 exhibits many features characteristic of (low ionization) Broad 
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Absorption Line (BAL) QSOs: high optical polarization (20% at 2800A, Smith et al. 1995), strong 
Fe II emission (extremely strong in Mrk 231) (Boksenberg et al. 1977), red optical colors, and a large 
IR/optical flux ratio. ASCA and ROSAT observations of Mrk 231, analyzed by Iwasawa (1999) and 
Turner (1999), revealed an extended, soft X-ray component, with a luminosity L ~ 10 ergs , 
and a hard component exhibiting considerable absorption, with a 2-10 keV luminosity of only 
Lk 2 x 10 42 ergs _1 . If this is the intrinsic X-ray luminosity, then an AGN is unlikely to dominate 
the energetics unless Mrk 231 is X-ray-quiet. However, the signal-to-noise of the ASCA spectrum 
was too low to determine whether the observed X-ray emission is largely reflected, in which case 
the true X-ray luminosity could be much higher. 

In this letter, we present an analysis of a much deeper ASCA observation of Mrk 231 than that 
presented by Iwasawa. We find strong evidence for a powerful, buried AGN in this ULIRG which 
we observe only through reflection and scattering. Indeed, we show that the AGN in Mrk231 may 
well be the dominant energy source. 

2. Observation and Data Reduction 

Mrk 231 was observed by ASCA on 1999-Nov-10. The SIS data were cleaned in order to remove 
the effects of hot and flickering pixels and subjected to the following data-selection criteria : i) the 
satellite should not be in the South Atlantic Anomaly (SAA), ii) the object should be at least 7 
degrees above the Earth's limb, hi) the object should be at least 25 degrees above the day-time 
Earth limb, and iv) the local geomagnetic cut-off rigidity (COR) should be greater than 6GeV/c. 
We also applied a standard grade selection on SIS events in order to further minimize particle 
background. The GIS data were cleaned to remove the particle background and subjected to the 
following data-selection criteria : i) the satellite should not be in the SAA, ii) the object should be 
at least 7 degrees above the Earth's limb and hi) the COR should be greater than 7GeV/c. After 
these data selection criteria were applied, there are 96.3 ksec of good data per SIS detector and 
99.2 ksec of good data per GIS detector. 

Full-band images were extracted from each of the four ASCA detectors. Two point-like X-ray 
sources were found in the SIS image separated by 2 arcmin, with the dominant source (containing 
about 75% of the counts) coincident with Mrk 231. The GIS has insufficient spatial resolution 
to distinguish these sources. We extracted the SIS spectrum of each source and found them to 
have very similar shapes. Given that the fainter source does not dominate at any energy, and that 
detailed spectroscopy on these small scales is difficult due to the energy-dependent point spread 
function of ASCA, we chose to extract and analyze the combined spectrum of these sources (which 
is dominated by Mrk 231). 

Source photons were extracted from a circular region centered on the two sources with radii of 
3 arcmin and 4 arcmin for the SIS and GIS respectively. Background spectra were obtained from 
blank regions of the same field (using the same chip in the case of the SIS). No temporal variability 
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was observed. In order to facilitate x 2 spectral fitting, all spectra were rebinned so as to contain 
at least 20 photons per spectral bin. In order to avoid poorly calibrated regions of the spectrum, 
the energy ranges considered were 1.0-10 keV for the SIS detectors, and 0.8-10 keV for the GIS 
detectors. Note that we use a lower-energy cutoff for the SIS that is considerably higher than the 
'standard' 0.6 keV cutoff in order to avoid the effects of "residual dark current distribution", or 
RDD, which is known to plague recent ASCA observations. 

3. Spectral Analysis 

Iwasawa (1999) finds evidence for soft X-ray emission from thermal plasma and measures a 
temperature of kT = 0.98lo'gg keV and an abundance of Z = 0.08+0 05 ^Q- Since we deem our 
soft (E < 1 keV) SIS data to be unusable, we cannot constrain the nature of the thermal plasma 
component. Therefore, in all that follows, we constrain the parameters of the thermal component 
(kT and Z) to the 90% error range found by Iwasawa. 

Guided by Iwasawa (1999), we initially fit our SIS & GIS spectrum (Figure 1) with a model 
consisting of this thermal component and a power-law tail with photon index T (representing the 
AGN emission), all absorbed by the Galactic column of Nu = 1-26 x 10 20 cm~ 2 . Although this 
spectral model produced an adequate goodness of fit (x 2 /dof = 304/305), the resulting photon 
index was extremely flat (r = 0.69 ± 0.12). In particular, it is much flatter than the intrinsic 
spectrum of any known AGN (which usually have T = 1.5 — 2.2). Adding intrinsic absorption to 
the power-law component increases the permitted photon index by Ar ~ 0.4 (with Ax 2 = 7 for 1 
new parameter) — insufficient to bring T into the acceptable AGN range. 

Since this AGN is likely to be highly obscured, we also examined classes of models in which the 
central power-law is completely absorbed by a Compton-thick screen, and is only seen via scattering 
and Compton reflection. In practice, we used the PEXRAV models within the xspec package, 
interpreting finite reflection fractions as representing mixed reflection/scattering (i.e., scattering 
of non-reflected as well as reflected continuum into our line of sight). A model consisting of this 
reflected/scattered component and the thermal plasma produced a good fit (x 2 /dof = 284/304) 
and implied a reflection-dominated spectrum (with the reflection fraction 1Z > 5800). However, the 
implied intrinsic power-law was much steeper (r = 3.01+q J; 2 ;) than the normal AGN range. 

Physically acceptable results are produced when absorption of the reflected/scattered com- 
ponent is included. While this absorption is not statistically required, including a column of 
Nu ~ 3 x 10 22 cm~ 2 extends the range of allowable photon indices into the typical AGN range. 
In order to better constrain the fits, we hereafter assume that the hard X-ray photon index of the 
source is typical of radio-quiet AGN, T = 1.8. We also include a narrow iron Ka emission line of 
cold iron which is expected to accompany the X-ray reflection. The best fitting parameters (with 
90% error range) are then: 1Z = 8. 7+54, N-r = 2.9t 2 'o x l^ 22 cm ~ 2 > an d an iron line equivalent 
width of Wxa = 290l^ eV (with x 2 /dof = 280/303); the resulting model spectrum is shown in 
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Figure 2. We have experimented with fits assuming partial covering of the source or ionized rather 
than neutral absorbers; neither the quality nor the nature of the best fits differ in any significant 
way from the above. This is to be expected since the hard X-ray spectrum largely decouples in the 
analysis from the soft, absorbed, spectrum. 

We also investigated the potential presence of an unabsorbed power-law component (corre- 
sponding to AGN emission that has been scattered at sufficiently large distances to be free of any 
substantial absorption). No such component was required by the data, and limits were set on the 
luminosity of any such component of < 1.1 x 10 42 ergs -1 . 



4. Discussion 

In his analysis of previous ASCA data, Iwasawa (1999) was unable to distinguish between a 
very flat, unabsorbed power-law and a heavily absorbed, typical (r = 1.8) AGN spectrum. With 
the considerably better S/N provided by our ASCA observation, we can show unequivocally that 
no amount of neutral absorption can steepen the intrinsic photon index into the normal AGN 
range. We are therefore led to conclude that there is a Compton-thick screen blocking the X-ray 
source along our line of sight, and we see the AGN only by reflection and scattering of the X-rays 
into our line of sight (Figure 3). Furthermore, unless the spectrum of the hard X-ray emission is 
unusually steep, the reflected/scattered emission is fairly heavily absorbed, by a neutral column 
N H « 3 x 10 22 cm" 2 . 

We can estimate the intrinsic luminosity of the X-ray source in two ways. When we subtract the 
contribution from the thermal emission, the observed 0.5 — 10 keV luminosity is 2.0 x 10 ergs -1 . 
Corrected for neutral absorption, this becomes 3.2 x 10 ergs -1 , and with the reflection fraction 
TZ set equal to zero, the scattered luminosity is L sc = 2.4 x 10 ergs -1 . (Note that this implies 
that about three-quarters of the observed X-ray luminosity is due to scattering of X-rays into our 
line of sight, rather than reflection, thus reducing the equivalent width of the Fe Ka line from the 
value of ~ 1 keV expected for a reflection-dominated spectrum.) Hence the intrinsic luminosity of 
the X-ray source is given by 

L intr = 2.4 x 10 43 (j^j ergs" 1 (1) 

where the scattering fraction / sc is expected to be on the order of 1% for electron scattering, 
assuming that the scattering column density is of the same order as the inferred absorption column 
density, i.e., a Thompson optical depth of order 0.01. 

Alternatively, we can estimate L mtr from the reflected fraction and modeling of the reflection 
process. From the above numbers, the observed luminosity in the reflection component is L re fl = 
8.0 x 10 41 ergs -1 . Using xspec, we conclude that the product of L- mtI and the reflection fraction 
Acfl is 

L intr / refl = 2.0 x 10 43 ergs -1 (2) 
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and so 

L intr = 2.0 x 10 44 ( ergs" 1 (3) 

V/refl/ 

where we have normalized to 0.1 under the assumption that only a small fraction of the reflecting 
surface is visible to us; for the Compton-thick Sy 2 NGC 1068, Iwasawa, Fabian Sz Matt (1997) 
infer / re fl ~ 0.05. These two estimates of the luminosity agree for a ratio / rG fl// sc ~ 10, i.e., for a 
reflected fraction / rc fl ~ 0.1 for a scattered fraction / sc ~ 1% which, as noted above, is the value 
expected if the scattering column is of the same order as the inferred absorbing column. (This ratio 
of /sc/ireA is not independent of the value 1Z ps 10 derived from the spectral fitting.) Thus we infer 
that /refi ~ 0.1, / sc ps 0.01, and L x ps 2 x 10 44 ergs _1 . This is 200 times larger than the value of 
Iwasawa (1999), and a factor ~ 20 larger than the luminosity estimated by Turner (1999); hence 
our analysis implies that the AGN is much more X-ray- luminous than previously thought. 

From our estimate of the intrinsic hard X-ray luminosity L x , we can estimate the contribution 
of the AGN to the bolometric luminosity of Mrk 231. The implied ratio of L x /Li, i ~ 2%. For 
comparison, the ratio of L x /Lb i for the sample of Sy 1 and radio-quiet QSOs of Elvis et al. (1994) 
for which X-ray data are available ranges from 1.95% to 14.5%, with the mean value decreasing 
somewhat with L bol , from 6.3% for L hol = 10 45 -10 46 (11 objects) to 4.3% for L hol = 10 46 - 2.2 x 10 46 
(7 objects); all the latter are classified as radio-quiet QSOs. Hence, within the uncertainties, the 
energetics of Mrk 231 are consistent with complete domination by the AGN; in any case, the 
fraction of the bolometric luminosity contributed by the AGN is almost certainly large ( £ 0.5). 
Extrapolating to high energies, the predicted 50 — 200 keV luminosity is -£50-200 ~ 2 x 10 44 erg s _1 , 
a factor of five below the Dermer et al. (1997) upper limit from OSSE. 

We also note that, with our estimate of the intrinsic X-ray luminosity, Mrk 231 is no longer 
an outlier in the significant correlations found by Lawrence et al. (1997) in the properties of AGN 
with extreme Fe II emission, in particular, those between the 2 keV - lfim spectral index ct\ x and 
the Fe II emission strength and the emission line FWHM. Since our estimate of the flux density at 
2 keV is ~ 270 times larger than the value used by Lawrence et al. (from Rigopoulou, Lawrence, 
& Rowan-Robinson 1996), we find a value of a- lx ps 1.2, smaller by 0.8 than the value inferred by 
Lawrence et al, which moves Mrk 231 into line with the bulk of strong Fe II emitters in their 
sample. Furthermore, the expected hard X-ray luminosity from the correlation between broad Ha 
and 2 — 10 keV emission (Ward et al. 1988), using the Ha flux measured by Smith et al. (1995, 
uncorrected for extinction), is L x <L 2 x 10 44 ergs _1 , in excellent agreement with the inferred X-ray 
luminosity. 



5. Summary 

A moderately long ASCA observation of the infrared- luminous galaxy Mrk 231 shows a very flat 
hard X-ray spectrum, which strongly suggests that the central AGN is hidden behind a Compton- 
thick screen, and is observed in X-rays only by scattering and reflection of the radiation into our 
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line of sight. The inferred X-ray luminosity in the 0.5 — 10 keV band is L x « 2 x 10 ergs 
for a scattering fraction / sc ~ 1%, the value expected from the absorbing column density inferred 
from the X-ray spectrum. This is much larger than previous estimates, and the resulting ratio 
of L x /Lb i is typical of luminous Sy 1 and radio-quiet QSOs, indicating that the AGN dominates 
the bolometric luminosity of the galaxy, and confirming the suggestion of Smith et al. (1995) that 
Mrk 231 is most properly classified as a low-ionization BAL QSO. The increase in the inferred 
X-ray luminosity over previous estimates also puts Mrk 231 in line with the correlations between 
2 keV - 1/um spectral index and both Fe II emission strength and emission line FWHM found by 
Lawrence et al. (1997) for AGN with extremely strong Fe II emission. 
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Fig. 1. — ASCA SIS and GIS spectrum of Mrk 231. The top panel shows the data for all four 
instruments together with the best fitting folded model (i.e. the model has been folded through 
the response of the telescope/detector system). The bottom panel shows the contribution of each 
spectral bin to the overall \ 2 of the best fit. 
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Fig. 2. — Best fitting spectral model for Mrk 231, plotted in terms of photon flux. The re- 
flected/scattered emission from the AGN component dominates at hard X-ray energies, but is 
truncated at 1-2 keV by heavy neutral absorption. The line-rich thermal plasma spectrum domi- 
nates at soft X-ray energies. The strong iron-L complex at ~ 1 keV is clearly visible. 
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Fig. 3. — Schematic illustration of our inferred geometry for the X-ray emission from Mrk 231, in 
which we observe only scattered and reflected emission and do not have an unobstructed view of 
the nucleus, and there is a distribution of absorbing gas along our line of sight to both the reflected 
and scattered components of the X-ray emission. 



